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ABSTRACT 

Narrowband observations have detected no Lya emission within a 70 pMpc'^ vol- 
ume centered on the z — 2.168 quasar PKS 0424-131. This is in contrast to 
surveys of Lya emitters in the field at similar redshifts and flux limits, which 
indicate that tens of sources should be visible within the same volume. The 
observed difference indicates that the quasar environment has a significant in- 
fluence on the observed density of Lya emitters. To quantify this effect we have 
constructed a semi-analytic model to simulate the effect of a luminous quasar on 
nearby Lya emitters. We find the null detection around PKS 0424-131 implies 
that the minimum isothermal temperature of Lya emitter host halos is greater 
than 3.4 x 10^ K (68% level), corresponding to a virial mass of ~ 1.2 x 10^^ Mq. 
This indicates that the intense UV emission of the quasar may be suppressing the 
star formation in nearby galaxies. Our study illustrates that low redshift quasar 
environments may serve as a surrogate for studying the radiative suppression of 
galaxy formation during the epoch of reionization. 

Key words: cosmology: theory - galaxies: clusters: general, intergalactic 
medium, quasars: general ~ ultraviolet: galaxies 



1 INTRODUCTION 

Current models of reionization are constrained by obser- 
vation to begin at z; ~ 10 ( Komatsu et al.||2011[ ) and to 
have been completed by z ~ 6 ( Fan et al.|2006 1. The pop- 
ular picture for this process assumes that isolated and in- 
ternally ionized ultraviolet (UV) sources carved out bub- 
bles of ionized hydrogen (H ll) in the neutral intergalactic 
medium (IGM) . These bubbles grew in size and increased 
in number as the cosmic star formation rate increased 
and more UV sources illuminated the IGM. Eventually 
these bubbles overlapped until they pervaded all of space, 
leaving the entire IGM ionized and thus ending reioniza- 



for gas accretion onto a dark matter halo ( Dijkstra et al 



tion (for a review on reionization see Barkana & Loeb 



20011. The higher IGM temperature in these ionized re- 



gions raised the isothermal virial temperature required 
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2004 1 greatly increasing the critical mass required to form 



galaxies. This process of raising the minimum halo mass 
for galaxy formation - known as 'Jeans-mass filtering' 
- is thought to have played a crucial role in the transi- 
tion to an ionized IGM ( Gnedin|2000 1 . Learning how this 
mechanism works is therefore vital to our understanding 
of reionization ( Iliev, Shapiro fc Raga|2005 l. 



Observing Jeans-mass filtering during the epoch of 
reionization directly is not possible at present, making it 
difficult to determine its role with respect to completing 
reionization within the observed timeline. To constrain 
difi'erent mechanisms of galaxy formation during reion- 
ization with current instruments therefore requires a sur- 
rogate environment that can be readily observed, such as 
the dense and ionized regions around quasars. By collect- 
ing statistics about the number densities and masses of 
galaxies within the highly ionized and clustered regions 
around quasars we can get an observational handle on the 
effects of a highly ionized environment on these galaxies. 
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The goal of this paper is to compare the observed num- 
ber density of galaxies with a theoretical model in order 
to highlight any environmental effects introduced by the 
highly ionized environment. 

[Francis fc Bland-Hawthorn| ( |2004| hereafter |FBH04 1 



imaged the region around the luminous z = 2.168 quasar 
PKS 0424-131, looking for fluorescent Lya emission from 
clouds of neutral hydrogen (Hi). They used the Tau- 



rus Tunable Filter ( jBland-Hawthorn fc Jones 1998 1 on 
the Anglo-Australian Telescope to probe a volume of 
70pMpc^ centered on the quasar with three narrow- 
band (7 A FWHM) images tuned to rest frame Alyq at 
z = 2.161, z = 2.168, and z = 2.175. This technique 
provides low resolution spectra across the entire field of 
view, allowing any source of Lya emission to be easily se- 
lected by looking for dropouts in the three redshift bands. 
Based on surveys done at similar redshifts and account- 
ing for galaxy clustering, [FBH04]expected to see between 
6 and 25 fiuorescing hydrogen clouds of varying sizes and 
> 10 internally ionized Lya emitting galaxies (LAEs). 
However, their observations found no hydrogen clouds, 
nor any LAEs, leading them to the tentative conclusion 
that quasar induced photo-evaporation was destroying 
the clouds and preventing or suppressing the formation 
of stars in nearby galaxies. 

In this paper we construct a semi-analytic model to 
interpret the observation in |FBH04| A null detection of 
LAEs in these regions either means that galaxies do not 
exist near to the quasar (i.e. were destroyed or never 
formed) or that they are not detectable (their emission 
is obscured by dense patches of the IGM or by interstel- 
lar dust) at the epoch of observation. Our semi-analytic 
model includes the effects of Lya transmission and galaxy 
clustering, and can be used to explore the implications 
of radiative feedback on galaxy formation. We tailor this 
model to the environment of a luminous quasar, and pop- 
ulate it with star forming galaxies following a density pre- 
scription that is calibrated against luminosity functions 
from wide-field LAE surveys ( |Quchi et al.|2008] hereafter 
O08}. 

We describe the model in ^ and present comparison 
with observations in i|3] along with a discussion on how 
consistent a null detection is with this model. In !j4] we 
summarize our findings, argue that there is a strong indi- 
cation of the suppression of low mass galaxies within in 
the volume around PKS 0424-131, and suggest that fur- 
ther surveys of LAEs near quasars can yield constraints 
on the process of Jeans-mass filtering during reioniza- 



tion. We assume the standard WMAP7 cosmology (Ko- 
matsu et al.poTT| , (n„, ^a, f^b, h, ag, n,) = (0.27, 0.73, 
0.046, 0.70, 0.81, 0.96). Throughout we adopt the con- 
vention of specifying a distance as physical or co-moving 
by prepending a 'p' or 'c' to the distance unit (i.e. pMpc 
and cMpc). 



2 MODELING OF Lya GALAXIES IN 
QUASAR ENVIRONMENTS 

This section outlines the model used in our analysis. The 
model is divided into two parts, i) a calculation of the ob- 
served Lya luminosity and number density for galaxies of 



a particular mass, and ii) how these relationships are af- 
fected by the exotic environment of a nearby quasar. The 
observed Lya luminosity is determined by the intrinsic 



Lya luminosity {\2.1\ multiplied by the fraction of this 



luminosity that is transmitted through the IGM (j ]2.2[ |. 
The model also yields a UV luminosity {{ 2.3 1 allowing 
us to fit both Lya and UV luminosity functions to de- 
termine the free parameters in the model (j |2.4[ ). Finally, 
the model also accounts for the impact of the quasar on 



the IGM near the galaxies {{ 2.5 1 



2.1 Galactic Lya Luminosity 

To determine the intrinsic Lya luminosity correspond- 
ing to a given halo mass we follow the semi-analytic 



model outlined in Dijkstra, Lidz fc Wyithe (2007 here- 



after D07I. This model first determines the galaxy's star 
formation rate {Mi,) then uses this to calculate the re- 
sultant ionizing luminosity (Qh) and hence a Lya lumi- 
nosity (I/Lya). These relationships depend on the unitless 
free parameters for duty cycle (lifetime as a fraction of 
the Hubble time; edc), star formation efficiency (/,), and 
ionizing escape fraction (/obc). 

For a halo with mass m at redshift z we define the 
star formation rate to be: 



fDC thub (z) 



M, 



where {^l-t,/Qr„m) is the mass of baryons in the galaxy 
and er,cttiub{z) is the total time over which the galaxy 
was forming stars. The star formation rate is used to 
obtain the ionizing photon luminositj^ (Qh) following 
Kennicuttl (119981): 



Qh = 9.26 X 10^^A/,s"\ 



(2) 



which assumes a Salpeter IMF with stellar masses rang- 
ing 0.1 — 100 Mq. Assuming that two out of three ioniz- 
ing photons that do not escape the galaxy are converted 
to Lya through case-B recombination ( Osterbrock|1989 1, 
the final equation for the Lya luminosity is 



I/Lyc = 0.68 hp Ui^ya (1 — /osc) Qh ergs s 



(3) 



where is Planck's constant, i^Lya is the frequency of a 
Lya photon, and f^^a is the fraction of ionizing photons 
that escape the galaxy without being absorbed. Com- 
bining these equations gives us an expression for Llyq 
that depends on the total halo mass and redshift, and is 
proportional to the model parameters edc, /*, and /oac, 
which we assume to be mass independent. 



^ Our model is insensitive to the exact relationship between 
Qh to M» as any changes to the assumed evolutionary syn- 
thesis model are absorbed into the free parameters and f^^c 
during the luminosity function fitting. Thus the high-redshift 
low-abundance {Z 



Schaerer 



0.05 Zq), model for Qh in 
l |2003[ l produces the same results as the solar abundance model 
used in |Kennicutt| l |1998^ , which for simplicity is used for the 
remainder of this paper. 
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2.2 Lya Transmission in the IGM 

Observations of Lya-emitting galaxies are subject to res- 
onant absorption of Lya flux from H I atoms in the IGM. 



Gunn & Peterson ( 1965 1 showed that for high-redshift 
Lya? sources, even a modest neutral fraction of 10~^ can 
significantly reduce the number of transmitted Lya pho- 
tons along our line of sight, as the photons blueward 



of A, 



redshift through Lya resonance. Therefore the 



transmission of Lya photons through the IGM is a critical 
component of a model for the Lya luminosity function. 

To determine the impact that the IGM has on the 
transmission of Lya photons we make extensive use of 
the model outlined in §3 of [D07| The IGM is modeled 
following the halo infall calculations of 'Barkana] ( |2004| 
and interpretation by |D07[ who provide a calculation of 
the density profile (pigm) and velocity field (vigm) of the 
IGM as a function of halo mass and distance from the 



halo [compare with Dijkstra, Lidz fc Wyithe||2007 
tion (4)]: 



equa- 



"'^"\p[l + exp(2-./5r.„) 



l^IGM(r, rvir) = 



r < lOrvir, 
r > lOrvi,, 

r < lOrvir, 
r > lOr^ir, 



(4) 



where r is the distance from the halo, rvir and Udrc are the 
virial radius and circular velocity of the halo, and H{z) 
is the Hubble parameter at redshift z. 

The density profile is combined with the assumed 
photoionizing flux from the galaxy and the external UV 
background (Fbg) to get a distance dependent neutral 
fraction (xigm). These values are combined to give the 
total opacity of photons as a function of their wavelength 
A [compare with Dijkstra, Lidz fc Wyithe|2007 equation 
(8)]: 



r(A) 



dr Pigm (r) Xigm (r) crLy„ [A, «igm (r)] , (5) 



where ctlyc. is the absorption cross-section for the rota- 
tionally broadened Lya line, evaluated at wavelength A, 
and blueshifted by the velocity field Vigm • This opacity is 
then convolved with an assumed IGM density fiuctuation 
distribution ( Miralda-Escude, Haehnelt fc Rees|2000 1 to 
account for any clumpy overdensities along the line of 
sight. The resulting function {e~^)(A) gives the fraction 
of photons emitted at wavelength A that are transmitted 
through the IGM without being scattered out of the line 
of sight. 

The quantity of interest is the total fraction of the 
Lya line transmitted through the IGM: 



JdA(e-)(A)J(A) 
/dAJ(A) ' 



(6) 



where J(A) is the flux of the galaxy as a function of wave- 
length. A galaxy with an intrinsic luminosity above the 
detection limit can be pushed below detectability for a 
sufficiently small value of TTgm. 

Fig.rn shows {e~'^)(A) across the width of the rest 
frame Lya line for a galaxy with mass IO^^Mq in mean 
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Figure 1. The fraction of doppler-shifted Lyo photons trans- 
mitted through the IGM from a galaxy with total mass 
10^^ Mq as a function of wavelength for field conditions at 
redshifts (from top to bottom) z = 2.2, 3.1, 3.7, and 5.7. The 
integrated transmission across the broadened Lyo for each red- 
shift is (from top to bottom) TIgm = 0.74, 0.48, 0.31, 0.13. For 
reference the vertical line is the intrinsic Alyq . 



conditions at four redshifts. The assumed shape of the 
broadened Lya line is Gaussian with its width set by 



DO' 



the assumed dark matter circular velocity following 
The four curves correspond to the redshifts (top to bot- 
tom) 2 — 2.2, 3.1, 3.7, and 5.7, with integrated transmis- 
sion fractions (top to bottom) Tigm = 0.74, 0.48, 0.31, 
and 0.13. For reference the vertical line is ALya,. The 
sharp drop in transmission around 1216 A is caused by 
the blueshifted Lya resonance as seen by photons escap- 
ing the galaxy through infalling hydrogen gas from the 
IGM. Between 1215 and 1216 A the galaxy's internal ion- 
ization and the ionizing background decrease the opac- 
ity of the nearby IGM to Lya photons causing the blue- 
ward rollup in (e~'^)(A), and as the ionizing background 
and UV mean free path increases with decreasing redshift 
this effect becomes more pronounced. Note that the total 
width of the broadened Lya line changes as the circular 
velocity of the galaxy changes with redshift, resulting in 
different wavelength spans. 

The shape of the transmission curve changes with 
redshift as the model components evolve, with the biggest 
contributions coming from the decreasing mean hydrogen 
density which evolves from 5.8 x 10~^ cm~^ at z = 5.7 to 
6.1 X 10"*^ cm"^ at = 2.2, and the external UV back- 
ground which evolves from 0.34 x 10~^'^ s~^ a.t z = 5.7 to 



1.20 X 10"^^ s"^ at 2 = 2.2 ( [Bolton fc Haehnelt||2007 l. 
The decrease in IGM density and increase in Fbg allows 
for the IGM transmission fraction to increase by a factor 
of six between redshifts 5.7 and 2.2. This illustrates the 
sensitivity of the Lya transmission to these two environ- 
mental factors, and motivates the idea that the quasar 
environment may be significantly reducing the transmis- 
sion of nearby Lya emitters. 

Fig. [2] shows (e~'^)(A) convolved with an assumed 
Lya line profile to simulate the continuum subtracted 



^ [Dijkstra fc Wyith"e| (jSoToJ show that galactic outflows of H I 
modify the Lya spectral line shape decreasing the impact of 
the IGM on transmission, implying that our simple Gaussian 
model is a conservative calculation of Tigm. 
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Figure 2. Lyo line shapes for a 10 Mq galaxy at 
in various environments. The dotted curve is the original Lya 
line shape and the grey dashed curve is the transmitted line 
shape assuming a mean IGM density and UV background. 
The solid black (red) curves show the impact of an increased 
IGM density and UV flux from the nearby quasar with a host 
halo size of 2 x IO^^Mq (IO^^Mq). The dash-dotted curves 
include the increased IGM density but without the quasar's 
contribution to the incident UV flux. 



spectrum for a galaxy with total mass 10 Mq at z = 
2.2. The dotted curve is the original Lya line shape, and 
the grey dashed curve is the transmitted line shape as- 
suming a mean IGM density and UV background. The 
sharp feature around 3852 A is the blueshifted Lya res- 
onance as seen by photons escaping the galaxy through 
infalling hydrogen gas from the IGM. In the absence of a 
strong UV flux the entire line blueward of this shifted res- 
onance is scattered out of the Hne of sight and the trans- 
mission is lower than the 0.5 predicted by |Gunn fc Peter^ 



(19651. With a large enough UV flux from internal 



(the galaxy) and external (UV background) sources, the 
neutral fraction of the IGM in the vicinity of the galaxy 
is lowered considerably, and as a result, the transmission 
of the galaxy's intrinsic Lya luminosity is increased. 



2.3 Galactic UV Luminosity 

In addition to Lya luminosity our model also calculates 
a UV magnitude for each galaxy. This allows us to use 
an additional dataset to help constrain euc, /*, and f^sc 
at the cost of an additional fitting parameter /dust to 
account for UV-specific|^dust extinction. We assume this 
to be luminosity independent and fit this as an additional 
free parameter in our model ( [Dayal &: Ferrara|2011 \. Our 
fitted values for /dust match those found by |Bouwens et al.| 



(20091 



We use the relationship between star formation rate 
and rest-frame UV continuum luminosity (I/a.uv) as cal- 



culated by Kennicutt ( 1998 1 



La,uv = 7.1 X lO'^^fdu^tM. ergss"' Hz" 



(7) 



^ We do not explicitly account for dust extinction in our treat- 
ment of Z/Lya as any pre-IGM effects of dust on Llyq is ab- 
sorbed into /esc while fitting the luminosity functions. 



Combining this expression with Eqn. ([T| gives us the re- 
lationship between UV luminosity and the free model pa- 
rameters. To compare with the UV luminosity function 
in |O08[ this is converted to an absolute AB magnitude: 



M. 



-2.51ogio (L;,,uv) + 51.6. 



(8) 



2.4 Fitting UV and Lya Luminosity Functions 

To constrain our Lya model parameters edc, and 
/cBc we fit the 2 = 3.1 Lya and UV differential lumi- 



nosity functions presented in O08 (their figures 16 and 
22 respectively) using the luminosity model presented 
above. For each luminosity bin in the observed Lya and 
UV luminosity functions we determine the mass required 
to generate that luminosity by inverting Eqn. ([3| and 
Eqn. ([8| respectively. This mass is used to obtain the 
halo number density for that luminosity bin using the 
extended Press-Schechter mass function, dn/dm ( jSheth^ 
Mo & Tormen 20011. The halo number density is then 



converted to a galaxy number density by assuming only 
a fraction cdc of the halos are occupied by star form- 
ing galaxies in the observed epoch and therefore emit- 
ting a detectable I/Lyo and I/^.dv- This yields the follow- 
ing equations for the differential Lya and UV luminosity 
functions: 



> (7lGM X Z/L 



nuv(Muv) = Cdc 



dn dm 
dm dlog^gL ' 
dn dm 
dm dMuv ' 



(9) 
(10) 



where nLya is the number density of Lya galaxies with 
luminosities between log-^^L and logj^L + dlogj^L, and 
nuv is the number density of UV galaxies with AB mag- 
nitudes between Muv and Muv + dMuv . 

Both the Lya and UV luminosity functions are fit 
simultaneously by comparing to the data points in |O08| 
and computing a for both Lya and UV. The combined 
value Xtat = xL= + Xuv is minimized by varying /. , /„sc 
(Lya), /dust (UV) for fixed values of edc- The parameter 
tuc is highly degenerate in our model and so the resultant 
best fits for a range of values are presented in Table[][f| 

Fig.[3] shows our model fits to the z = 3.1 Lya and 
UV luminosity functions presented in |O08[ The fits cor- 
respond to the upper and lower bounds of edc = 0.4 
and 1.0 that are in agreement with the data. There are 
14 binned data points in total, 8 for the Lya luminosity 
function and 6 for the UV. This is fit with the 4 fitting 
parameters for a total of 10 degrees of freedom. We found 
that edc = 1.0 yields the best fit for the luminosity func- 
tions, with a reduced xLt of 0.82. This corresponds to 



* Our model assumes long-term, continuous, galactic star 
formation for simplicity. Recent work by [Sharp 



Bland-I 



|Hawthorn| | |2010[ | showed that star formation is impulsive on 
timescales less than < lO'' years and that for active galactic 
nuclei the timescales are much longer (> 10^ yr) corresponding 
to a > 3 cMpc light crossing time between cycles. This indi- 
cates that more sophisticated simulations will need a galactic 
duty cycle at least several times 10^ yr to be consistent with 
observation and that our modeling of a quasar's ionizing field 
as a continuous event is justified in our simulated volume size. 
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Figure 3. Lya (left panel) and UV (right panel) fits to z = 3.1 differential luminosity functions for two values of edc. The points 
and error bars are fromjOOS] figures 16 (left panel) and 22 (right panel). The best and worst fit values of cdc are plotted, edc = 1.0 
(solid curves) and edc = 0.4 (dash-dotted curves), with a reduced Xtot °f 0.82 and 1.90 respectively. Intermediate values of enc 
produce fits that lie between these two curves. 



Table 1. Best fit parameters for model fitting Lya and UV 
luminosity functions in |O08| 







/c=oxro-i 


/dust 




X ^ Xtot 


0.4 


7.7 


1.45 


0.22 


1.90 


4% 


0.6 


8.4 


1.15 


0.24 


1.32 


21% 


0.8 


8.9 


0.81 


0.26 


1.00 


44% 


1.0 


9.3 


0.56 


0.28 


0.82 


61% 



" d = 10 is the degrees of freedom in the fit. 

^ Probability that a random will be greater than Xtot ■ 

a probability of 61% that a randomized set of parame- 
ters would produced a larger Xtot, and we use this value of 
tuc for all calculations on the impact of a quasar environ- 
ment in the following section. The quality of fit decreased 
monotonically with edc and we found that below a 0.4 
the fit was in significant disagreement (< 4%) with the 
data. We will show in S|3]the exact choice of edc does not 
change our main result. 

As a further check we compared our calculated Lya 
equivalent widths (EW) to the mean observed at z — 3.1 
in |O08| F rom their spectroscopic and photometric sam- 
ples 



O08 



calculated a mean EW of 102 and 130 A respec- 
tively. Our best fit model has a mean EW of 120 A, m 
good agreement with these observations. Moreover, be- 
cause EW oc 7Igm(1 — /Gac)//duat this is independent ev- 
idence of a robust fit of our parameters to observation. 
Our range of /duut = 0.22 to 0.28, also matches nicely 
with recent results in [Blanc et al.| ( |2010[ ) who found this 
parameter to be between 0.20 to 0.30. 

We note that all of the constraints on our parameters 
are obtained by fitting to the z = 3.1 luminosity functions 
observed in |O08| but are then applied to an environment 
at z = 2.2. Current Lya luminosity functions at 2 < 3 
still suffer from cosmic variance (see [Blanc et al.||2010} 
and references therein) and there exist no accompanying 
UV luminosity functions at these lower redshifts. There- 
fore fitting to data from z = 3.1 is the best that can be 
done to constrain the free parameters in the model at 
lower redshift. 

It is therefore important to note that recent observa- 
tions do show evolving LAE properties with redshift, in 



particular the EW distribution at z = 2.3 (Nilsson et al. 
,2009. ) is described by a steeper exponential function than 
at z = 3.1 ( [Gronwall et al.[|2007| ) yielding a base EW at 
z — 2.3 that is roughly two-thirds of the value found at 
2 = 3.1. This suggests that Lya photons are more diffi- 
cult to detect relative to the UV continuum from z — 2.3 
galaxies than those emitted from z = 3.1. Therefore, to 
match the smaller 2 = 2.3 mean EW would require either 
an increase in our /dust and/or /esc parameters. Without 
a matching set of Lya and UV luminosity functions at 
this lower redshift we can only conjecture on exactly how 
this would affect our free parameters. Changing only /osc 
to match this reduced EW would have the largest effect 
on our results. We therefore ran our simulations with this 
lower /esc and found that the detected number of LAE 
was lowered by 35%, but that this did not change our fi- 
nal results. For the rest of the paper the fitted parameters 
in Table[T]are used. 



2.5 Quasar Influenced Transmission 

The additional component of our model beyond the work 
of |D07| concerns the environmental effects contributed by 
a nearby quasar, including an enhanced piom and Fbg- 
These have a significant impact on the Lya line shape, 
and are computed semi-analytically from the quasar's ob- 
served B-band apparent magnitude, ttib, and its redshift, 
2q. To determine both the increased ionization rate and 
IGM density in the vicinity of the quasar we first calcu- 
late its B-band luminosity Lb . Assuming that the quasar 
is shining at the Eddington limit we then find the black 
hole mass required to emit at the observed luminosity, 
and in turn calculate the expected host halo size using 
the conversion in [Wyithe fc L oeb| ( |2005t . 

To calculate the effective density of the IGM in the 
vicinity of a galaxy, we first write Eqn. (|4| as a density 
excess relative to the mean density p: 



'^p{r, T-vir) = Pigm(t-, rviO - p. 



(11) 



This density excess is then added to the underlying den- 
sity contribution of the quasar to get the total effective 
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density profile used in calculating TIgm- For a galaxy at 
a distance Tq from a quasar with viral radius rQ ^i, , the 
combined density of the local IGM a distance r from a 
galaxy with virial radius rvi, is: 



Pgai(r-,rvir) = piGM(rQ,rQ,vir) + Ap(r,rvir). 



(12) 



With this formulation, a galaxy just outside the virial 
radius of the quasar (rq = rq ^ir) has an effective den- 
sity that ranges from a maximum value of 39 p at the 
galaxy's virial radius (r = r^ii), to a minimum of 20 p at 
large distances from the galaxy (r > 30rvir). For a galaxy 
well away from the quasar (rg > SOrq^vir), the effective 
density ranges from 20 p at the galaxy's virial radius to 
the mean IGM density p at large distances. 

To calculate the ionization rate of the quasar we fol- 



low the prescription of Schirber & Bullock (20031 and 
determine the quasar's flux at the Lyman limit (Jll) 
given its observed Lb. This, along with the assumed UV- 
continuum slope of the quasar and the cross-section of 
hydrogen, is used to calculate the number of hydrogen 
ionizations per second (Fq) as a function of quasar lumi- 
nosity and distance from the quasar: 



rQ(r,LB) = 



12.0 / Jll (r, Lb) 

X exp(-r/r„fp) 10"^%"\ (13) 



where auv is the UV-continuum slope of the quasar, and 
r„,fp is the inferred mean free path of UV photons attenu- 
ating the ionizing flux ( Faucher-Giguere et al.|2008| ). We 
generate the ionization field for each galaxy using a two 
component model consisting of the distance dependent 
quasar flux in Eqn. ( |13[ ) and a constant ionizing back- 
ground]^ 

As an example. Fig. [2] shows the environmental ef- 
fects on transmission for a 10^^ Mq galaxy 500pkpc 
from a nearby quasar. The black curves assume a central 
quasar matching the properties observed in|FBH04] with 
an apparent B-band magnitude of 17.6, a correspond- 
ing B-band luminosity of 4.4 x 10^'^ Lq g, and inferred 
black hole and host halo masses of 7.6 x 10^ Mq and 
2 x 10^^ Mq respectively. The UV flux generated by the 
quasar in the vicinity of the galaxy very nearly removes 
the IGM's effect on the Lya line, yielding a transmission 
of 0.94, seen in the solid black curve. The black dash- 
dotted curve shows how the enhanced density of the IGM 
decreases the transmission when the UV flux is not en- 
hanced by the quasar, resulting in a final transmission of 
0.33. This is compared to the dashed curve showing the 
Lya line assuming a mean IGM and background UV fiux, 
with a total transmission of 0.74. 

For comparison, the red curves show the same galaxy 
in the vicinity of a much smaller and less luminous 



^ This assumption was compared with a numerical model that 
treated each galaxy discretely within the quasar's zone of influ- 
ence (rq < SOrq^vir) and used a smoothed ionizing luminosity 
from 30rq_vir to 10r„tp. The sum of the galaxy contribution 
and the smoothed background was flat across the entire vol- 
ume, with the clustering of galaxies pushing the flux up by 
~ 2% at the center of the simulation volume. This justifies 
the simple two component model of the effective ionizing flux. 




0.3 1.0 
Radius [pMpc] 



Figure 4. Averaged Lya transmission as a function of radius 
from a central quasar with host halo size (black curves) 2 X 
lO'^^ Mq and (red curves) 10^'^ M0- The grey dashed curve is 
the average transmission assuming a mean IGM density and 
UV flux. The solid curves include the density and UV flux 
from the nearby quasar and the dash-dotted curves include 
the increased density but with no UV boost. 



quasar, with an apparent _B-band magnitude of 23.0, a 
corresponding _B-band luminosity of 3 x 10^^ Lq g, and 
inferred black hole and host halo masses of 5.3 X IO'^Mq 
and 10^^ Mq respectively. For this smaller quasar the UV 
fiux is only sufficient to compensate for the increased den- 
sity of the IGM, as seen in the solid red curve which 
roughly matches the dashed background curve with a to- 
tal transmission of 0.71. Because the quasar has a smaller 
host halo the density boost to the IGM is not as strong at 
the same radius, which allows for a transmission of 0.56 
in the absence of the quasar's UV fiux seen in the red 
dash-dotted curve. 

The exact shape of the radiation fields emitted from 
a quasar is an open and thorny question. If the emissions 
are powerful and tightly coUimated the ionization is still 
likely to be diffused through the volume in some way, 
either by scattering off of dust grains or electrons via 
Thompson scattering ([Bland-Hawthorn, Sokolowski &:| 
CecIIlpgT] [Sokolowski, Bland-Hawthorn fc Cecn||1991| ). 
Moreover, if the beam moves with time then a time- 
averaged UV field must be considered. The impact of 
these details can be averaged over by analyzing large sur- 
veys containing observations like that in FBH04, This al- 
lows for the limiting cases to be described by a dilution 
factor /dii representing the fractional solid angle that the 
radiation field emitted by the quasar permeates, where 
/dii ~ 1 corresponds to a perfectly isotropic quasar and 
/dii = to a quasar emitting no ionizing radiation. For 
the purposes of this paper, in which we analyze the re- 
sults of a single observation, we will only look at these 
two extremes. 

Fig.|4] shows the Lya transmission averaged over ob- 
servable galaxies, assuming thcJFBH04jdetection limit, as 
a function of distance from the quasar a.t z — 2.2. As in 
Fig.[2] the black and red lines represent a central quasar 
with host halo mass 2 X 10" Mq and 10^^ Mq respec- 
tively, the solid and dash-dotted lines show the effects of 
the quasar with {fan — 1) and without (/dn = 0) the en- 
hanced ionizing flux. The dashed grey line is the average 
transmission for a mean IGM. For each quasar the virial 
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radius was used as the minimum radial distance in the 
model to avoid confusion in halo occupation; the larger 
and smaller quasars have a minimum radius of 280 pkpc 
and 100 pkpc respectively. 

For the smaller quasar we see that Wigm and Fq com- 
pensate for one another, leading the a radial transmission 
curve that is nearly equal to that of mean IGM. The 
larger quasar has a much larger UV flux and is able to 
ionize the denser IGM efficiently, leading to a transmis- 
sion roughly 20% higher than the mean. Assuming that 
a quasar's ionizing flux is isotropic, this implies that the 
transmission in the vicinity of a quasar at z = 2.2 is com- 
parable to the mean IGM for modest sized quasars and is 
higher for the very luminous. Therefore if TIgm is the pri- 
mary factor influencing detection of Lya emission from 
galaxies near to a luminous quasar, we then would ex- 
pect to see an increase in the number of observable Lya 
galaxies in the vicinity of the quasar observed by |FBH04l 



3 COMPARISON WITH OBSERVATIONS 
AROUND PKS 0424-131 

Having described the model, we are now in a position 



to interpret the observations of PKS 0424-131 (FBH04I 



We construct a physical volume around PKS 0424-131 
matching the observed volume, and fill it with galaxies 
following our semi-analytic prescription. We do this by 
first building a series of concentric spherical shells cen- 
tered on our quasar's modeled halo, distributed evenly 
in log-space. The shell radii range from the quasar halo's 
virial radius out to a radius that encloses the observed 
volume. These shells are then filled with dark matter 
hales over a range of masses as given by the Press & 



Schechter ( 1974 1 mass function (as modified by Sheth,^ 
Mo fc Tormen||2001 1, enhanced by the dark matter two- 
point correlation function with the linear bias of |Mo fc] 
White (19961. We set the minimum halo mass of the 
model to produce an intrinsic luminosity I/Lya at the 
survey detection limit. These halos are populated with 
galajcies as in §2.4[ by assuming all the halos are popu- 
lated but that only a fraction eDc of these galaxies are 
star forming (and thus detectable) in the observed epoch. 
We use the best fit parameters based on our analysis of 
the field luminosity functions. 

Thus the average number of halos Afg.^i{m,r), with 
mass in the range m±dm/2, with a separation of r±dr/2 
from a central quasar of mass M at redshift z is: 



A/'gai(m,r) = eDc(47rr^dr) 



dn(m, z) 
dm 



[l + ar,m,M,z)], (14) 



where dn/dm is the halo mass function, and ^ is the 
linearly biased two-point correlation function. This is 
the prescription found in §2.1.1 of Dijkstra et al.| ( [2008 l 
[Eqn. ( |14[ ) is a modification of their equation (1)], which 
describes the clustering of galaxies around a central dark 
matter halo but assumes a non-linear treatment of the 
separation bias. As we are probing much larger separa- 
tions a linear treatment is sufficient for our purposes. 
|FBH04| surveyed a field of view measuring 4.2pMpc 



by 2.0pMpc over a redshift range of Az = 0.015 cen- 
tered on the z = 2.168 luminous quasar PKS 0424-131. 
This range corresponds to a physical depth of 6.5pMpc 
at this redshift and an effective depth of 8.5 pMpc when 
accounting for peculiar velocity induced redshift distor- 
tion ( Kaiser|1987 L This yields an effective survey volume 
of 70pMpP^ 

Their detection limit for Lya emitting galaxies 
within this volume was 9.6 x 10~^* ergscm"^ s~^, corre- 
sponding to an ideal (Tigm = 1) minimum mass between 
3.8 X 10^° Mq (enc = 0.4) and 7.1 x 10^'' Mq (edc = 1-0). 
In our simulation we consider galaxies with host halos 
ranging from this minimum detectable mass to a maxi- 
mum of 10^^ Mq for completeness, and use 500 mass bins 
distributed equally in log-space to span this range. Using 
the inferred mass of 2 x lO" Mq for PKS 0424-131 de- 
scribed in §2.5| we calculated a virial radius of 280 pkpc 
which we used for the innermost spherical shell. The out- 
ermost shell was given a radius 4.84 pMpc, fully enclos- 
ing the rectangular observed volume. We use 100 radial 
shells to span this range of radii, distributed equally in 
log-space. 

For each of these radial shells we calculate the 



quasar's UV flux from Eqn. ( 13 1 and combine it with the 



external UV background to give an effective ionizing rate 
for the shell, and determine the quasar-enhanced IGM 
density and IGM neutral fraction detailed in §2.5| For 
each of the mass bins within each shell we determine 
the transmission given the quasar-influenced conditions 
and determine which masses have a detectable luminos- 
ity post-transmission. For the masses that are above this 
threshold we sum the average number of galaxies ex- 



pected at this separation given by Eqn. (141. The total 



number of galaxies for each radial bin are added to give 
a final value for expected average enclosed galaxies within 
the survey volume. 

Because our model is generated from spherical shells 
and the effective observed volume around PKS 0424-131 
is rectangular with dimensional ratios of roughly 4:2:1, 
we determine for each shell the fractional volume that 
is located outside the observed rectangular volume and 
reduce that shell's contribution to the enclosed galaxy 
count accordingly. This properly treats the radially de- 
pendent clustering and UV flux in both the long and short 
axes of the observed volume. 

Fig. [5] shows the results for this process over the al- 
lowed range of edc- The solid black curve shows that the 
average number of galaxies that should be visible within 
the volume surveyed by |FBH04| given the fiducial model 
including the quasar is between 70 and 97. Even if the 
ionizing flux from the quasar itself is omitted, as seen 
in the dash-dotted black curve, the average number of 
visible galaxies is still between 49 and 67. This can be 
contrasted with the grey dashed line which gives the ex- 
pected number of galaxies assuming a mean density and 
IGM with no quasar, which ranges between 21 and 26, 
and the long-dashed black line which gives the galaxy 
count assuming a constant Tigm = 0.5, ranging between 
40 and 54. The clustering of galaxies around the quasar, 
combined with the higher than background transmission 
rates seen in Fig.[4j pushes the expected number of galax- 
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Figure 5. Predicted number of galaxies contained within 
the Francis-Bland-Hawthorn volume as a function of mod- 
eled duty cycle for a quasar with host halo size (black curves) 
2x lO^'^ Mq and (red curves) 10^-^ Mq. The grey dashed curve 
is the number of galaxies expected in the field. The long- 
dashed lines show the enclosed count assuming a constant 
TiGM = 0.5. The solid (dash-dotted) curves include cluster- 
ing and transmission effects including (excluding) the boosted 
UV flux from the quasar. The upper and lower dotted curves 
for each quasar show the enclosed count excluding galaxies 
with T^i, < 5 X lO^K (< 6.6 x IQIOMq) and T^i, < 10** K 
(< 1.9 X 10^^ Mq) respectively. 



ies up well above background levels, in contrast with ob- 
servation. 

We have also repeated the above calculation for the 
same observable volume and redshift but for a much less 
luminous central quasar. The mass of the dark matter 
halo hosting the quasar was set to 10^^ Mq and the inner- 
most spherical shell's radius to this quasar's lower virial 
radius of 100 pkpc, keeping the rest of the conditions the 
same. In this case the lower IGM density and lower UV 
flux conspire to give close to the same number of galaxies 
for conditions when the quasar UV flux is considered and 
when it is neglected. The solid and dash-dotted red lines 
show the expected galaxy count in these cases, which 
varies between 34 and 49 over the allowed range in edc- 
The long-dashed red line gives the galaxy count assuming 
a constant 7igm ~ 0.5, ranging between 25 and 33. 

Thus our model predicts a number of galaxies that is 
far in excess of the null-detection in |FBHd4| This points 
to either an additional suppression of the Lya signal that 
has not been modeled, or else to some mechanism for the 
suppression of star formation. To investigate this possi- 
bility we introduce a simple cutoff in mass below which 
galaxies are unobservable. The dotted curves in Fig. [5] 
show the results from these imposed mass cuts, which 
correspond to T^i, < 5.5 x lO'^ K [< 7.6 x 10^° Mq] and 
Tvi, < 8.4 X lO'' K [< 1.4 X 10" Mq]. In our model these 
halos yield the Lya flux limit with a constant Tigm = 0.5 
for cdc = 0.4 and Edc = 1-0 respectively. The top black 
(red) dotted curve shows the enclosed count excluding 
galaxies with T^i, < 5 x lO'^ K [< 6.6 x 10^° Mq] and the 
bottom black (red) dotted curve shows the enclosed count 
excluding galaxies with Tvi, < lO'' K [< 1.9 x 10"Mq]. 
The range of enclosed galaxies for the Tvi, < 5 x 10^ K cut 
runs from 46 to 97 for the larger quasar (upper black dot- 
ted curve) and from 28 to 49 for the smaller (upper red 
dotted curve). For the Tvi, < 10^ K cut the enclosed count 



runs from 16 to 41 for the larger quasar (lower black dot- 
ted curve) and from 9 to 24 for the smaller quasar (lower 
red dotted curve). 

The mean number of galaxies required for the |FBH04| 
null detection around PKS 0424-131 to be consistent with 
our model at the 68% level is ~ 2.3, and at the 90% 
level is ~ 1.1. Thus, for edc = 0.4 (edc = 1.0) our model 
must exclude masses below 1.2 X 10^^ Mq (2.5 x 10^^ Mq) 
in order to bring the mean number of galaxies down to 
a value consistent with observations at the 68% level, 
and must exclude masses below 2.1 X 10^2 Mq (4.2 x 
10^^ Mq) to be consistent at the 90% level. This means 
that the most conservative virial temperature consistent 
with observation is T,;, > 3.4 x lO'' K (~ 1.2 x 10^^ Mq). 
This result may imply considerable radiative suppression 
of galaxy formation by the nearby quasar. 



3.1 Additional Comparisons 

Several other studies have investigated the population of 
Lya emission near luminous quasars. |Cantalupo, Lilly fc| 
Porciani (20071 reported a detection of 13 Lya sources 
clustered around a quasar at z = 3.1 in a ~ 25pMpc"^ 
volume (sparsely sampled from a larger ~ 200pMpc'^ 
volume), two of which they suggested were hydrogen 



clouds fluoresced by the quasar's ionizing radiation. Our 
model, accounting for the sparse sampling, predicts ap- 
proximately 20 galaxies should be detected in this vol- 
ume. The possible externally fluoresced sources suggest 
there is enough ionizing flux around these high-redshift 
luminous to strongly impact neighboring Hi. 

[Kashikawa et al^ ( |2007[ ) conducted a deep wide field 
narrowband survey for Lyman break galaxies and LAEs 
around QSO SDSS J0211-0009 at z = 4.9. They sur- 
veyed ~ 830pMpc^ in the vicinity of the quasar and 
detected 221 LAEs. Our simulation predicted ~ 260 de- 
tectable LAEs for a similar volume, quasar, and redshift. 
[Kashikawa et ah] ( |2007[ ) found that while the observed 
Lyman break galaxies formed a distributed filamentary 
structure which included the quasar, the LAEs were pref- 
erentially clustered around the quasar while avoiding a 
vicinity of 4.5 cMpc from the quasar. This region was cal- 
culated to have a UV radiation field roughly 100 times 
that of the background, and this was posited to be sup- 
pressing the formation of LAEs. This spacial distribu- 
tion of LAEs reinforces the idea that the environment in 
the smaller volumes probed by our model - in the direct 
vicinity of the quasar - is where the majority of galaxy 
suppression takes place. 

Recent work by [Laursen fc Sommer-Larsen| ( |2011[ ) 
calculated TIgm using a sophisticated hydrodynamics 
code that accounts for radiative transfer and interstel- 
lar resonant scattering. Our results are consistent with 
the lower end of their Tigm values for corresponding red- 
shifts, implying our semi-analytic values are conservative, 
and thereby reinforcing our results. 



4 CONCLUSIONS 

In this paper we have presented a semi-analytic model 
that predicts the number of visible Lya emitting galaxies 
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around a central quasar, taking into account the quasar's 
impact on the local IGM density and neutral fraction. 
The free parameters of the model are determined by fit- 
ting Lya and UV luminosity functions taken from large 
field Lya surveys conducted by |Q08[ We use this model to 
interpret observations of a 70 pMpc^ volume centered on 
the z = 2.168 luminous quasar PKS 0424-131, in which 
no LyQ emission was detected ( |FBH04| . 

We find that this null detection of Lya emitting 
galaxies can only be explained in a scenario in which 
we introduce a simple cutoff in mass below which galax- 
ies are unobservable. In order for our model to be con- 
sistent with observations at the 68% (90%) level we 
need to exclude all masses below at least 1.2 x 10^^ Mq 
(2.5 X IO'^^Mq), corresponding to a virial temperature 
greater than 3.4 x 10^ K. This result may imply consid- 
erable radiative suppression of galaxy formation by the 
nearby quasar and motivates further observations of Lya 
emitters in the vicinity of luminous quasars. Understand- 
ing this process in more detail will ultimately help to con- 
strain the extent to which radiative suppression of galaxy 
formation took place during the epoch of reionization. 

Tunable filters (TP) are a powerful approach to prob- 
ing emission-line objects at any redshift ( Jones & Bland- 
Hawthorn|2001 1. In the coming decade, there are several 
facilities that are ideally suited to searches for emission- 
line objects, including the Osiris TF on the Grantecan 



10.2m (|Cepa et al.[20 03), MMTF on Magellan ( [Veilleux 
et al.|[20T0 l, and TF instruments under development for 
the NTT 3.5m and the SOAR 4m telescopes ( [Marcehn 
et al.||2008[ [Taylor et aL]|2010[ ). All of these instruments 
are well adapted to studying the impact of QSOs on their 
environs. We envisage that IR tunable filters operating 
with adaptive optics will be able to push to even higher 
redshifts and down to lower galaxy masses. 

While our simulations looked at mass cuts across 
the entire volume surveyed, these upcoming observations 
will provide the statistics required for spatial mapping 
of LAEs around their central quasar allowing for future 
models to constrain the critical ionizing flux required 
to disrupt galaxy formation. This is a fundamental un- 
known required by hydrodynamical N-body simulations 
of galaxy formation in the vicinity of AGNs and during 
the reionization epoch, and can be constrained using cur- 
rent generation instruments. 
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